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Pseudospin symmetry is a relativistic symmetry of the Dirac Hamiltonian with scalar and vector 
mean fields equal and opposite in sign. This symmetry imposes constraints on the Dirac eigen- 
functions. We examine extensively the Dirac eigenfunctions of realistic relativistic mean field calcu- 
lations of deformed nuclei to determine if these eigenfunctions satisfy these pseudospin symmetry 
constraints. 
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I. INTRODUCTION 


Pseudospin doublets were introduced more than thirty years ago into nuclear physics to accommodate an observed 
near degeneracy of certain normal parity shell model orbitals with non-relativistic quantum numbers (ny, £, 7 = €+1/2) 
and (n, — 1,4 + 2, j = L+ 3/2) where n,, £, and j are the single-nucleon radial, orbital, and total angular momentum 
quantum numbers, respectively [1, 2]. The doublet structure is expressed in terms of a “pseudo” orbital angular 
momentum, which is an average of the orbital angular momentum of the two orbits in the doublet, l = £ + 1, coupled 
to a “pseudo” spin, § = 1/2 with j = £ £ $. For example, the shell model orbitals (n,81/2, (nr — 1)d3/2) will have 


f= 1, (n;paj2, (nr — 1) fs/2) will have f= 2, for the two states in the doublet. Then the single-particle energy is 
approximately independent of the orientation of the pseudospin leading to an approximate pseudospin symmetry. 
These doublets persist for deformed nuclei as well [3]. The axially-symmetric deformed single-particle orbits with 
non-relativistic asymptotic quantum numbers [N,n3, AJQ = A+ 1/2 and [N, n3, A’ = A + 2]0' = A + 3/2 are quasi- 
degenerate. Here N is the total harmonic oscillator quantum number, n3 is the number of quanta for oscillations 
along the symmetry axis, taken to be in the z-direction, A and Q are respectively the components of the orbital and 
total angular momentum projected along the symmetry axis [4]. In this case, the doublet structure is expressed in 
terms of a ^pseudo" orbital angular momentum projection, A=A+ 1, which is added to a “pseudo” spin projection, 
jt = +1/2 to yield the above mentioned doublet of states with Q = A — 1/2 and € = A + 1/2. This approximate 
pseudospin “symmetry” has been used to explain features of deformed nuclei, including superdeformation [5] and 
identical bands [6, 7, 8, 9] as well. 

Although there have been attempts to understand the origin of this “symmetry” [10, 11], only recently has it 
been shown to arise from a relativistic symmetry of the Dirac Hamiltonian [12, 13] which we review in Section II. 
This relativistic symmetry implies conditions on the Dirac eigenfunctions [14] which we discuss in Sections II and 
III. These relationships have been studied extensively [14, 15, 16, 17, 18, 19] for spherical nuclei. For deformed 
nuclei, the relationships have been studied only in a limited way and primarily for the lower components of the Dirac 
eigenfunctions [20, 21, 22]. In this paper we shall test thoroughly these relationships between the upper and lower 
components of the two states in the doublet for realistic deformed relativistic eigenfunctions [23, 24]. 
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II. THE DIRAC HAMILTONIAN AND PSEUDOSPIN SYMMETRY 


The Dirac Hamiltonian, H, with an external scalar, Vs(r), and vector, Vy (7), potentials is given by: 
H 2 à: p B(M  Vs(P)) - WF) , (1) 


where &, B are the usual Dirac matrices, M is the nucleon mass, and we have set A = c = 1. The Dirac Hamiltonian 
is invariant under a SU(2) algebra for two limits: Vs(r) = Vy (r) + C; and Vs(r) = —Vv(r) + Cps where C,, C5, are 
constants [25]. The former limit has application to the spectrum of mesons for which the spin-orbit splitting is small 
[26] and for the spectrum of an antinucleon in the mean-field of nucleons [27, 28]. The latter limit leads to pseudospin 
symmetry in nuclei [12]. This symmetry occurs independent of the shape of the nucleus: spherical, axial deformed, 
or triaxial. 


A. Pseudospin Symmetry Generators 


The generators for the pseudospin SU(2) algebra, S; (i = x,y,z), which commute with the Dirac Hamiltonian, 
[Hs , Si] = 0, for the pseudospin symmetry limit Vs(r) = —Vy (T) + Cys, are given by [13] 


e n Si 0 E U, s; Up 0 
a= (4 er 0 4 (2) 


go-:p 
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where s; = c;/2 are the usual spin generators, o; the Pauli matrices, and U, = is the momentum-helicity 


unitary operator introduced in [11]. Thus the operators $; generate an SU(2) invariant symmetry of H,,. Therefore, 
each eigenstate of the Dirac Hamiltonian has a partner with the same energy, 


where k are the other quantum numbers and ji = +4 is the eigenvalue of Bas 
G DPS (A — 7 MPs (F 
5, BP (F) = pam (P) . (4) 


The eigenstates in the doublet will be connected by the generators S,—S,-c iSy, 


The fact that Dirac eigenfunctions belong to the spinor representation of the pseudospin SU(2), as given in Eqs. (4)- 
(5), leads to the conditions on the corresponding Dirac amplitudes that are explored in this paper and developed in 
the next Subsection. 


B. Dirac Eigenfunctions and Pseudospin Symmetry 


An eigenstate er (7) of the Dirac Hamiltonian H,,, Eq. (3), is a four-dimensional vector, 


PELO | apt a |: (6) 


where I; C) are the “upper Dirac components” and f; 9) are the “lower Dirac components". The superscript 
+ (—) indicates spin up (spin down). 


The connections between the Dirac eigenstates of the doublet (A = + 
relationships between the Dirac amplitudes in Eq. (6) [14], 


i) resulting from Eqs. (4)-(5) lead to 


jel ed, (Ta) 
RS = KLORA, (Tb) 
of) = Re =o), (c) 
and to first order differential equations 
ðo .a "TUN 0.0 "DT 

(um) = (rm) (8a) 
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'Thus pseudospin symmetry reduces the eight amplitudes for the two states in the doublet to four amplitudes. In the 
next Section we shall discuss these relations in Eqs. (7)-(8) for axially deformed nuclei. 


III. PSEUDOSPIN SYMMETRY FOR AXIALLY DEFORMED NUCLEI 


If the potentials are axially symmetric, that is, independent of the azimuthal angle $, Vs,v(r) = Vsv(p,z), p = 
va? + y?, then the Dirac Hamiltonian has an additional U(1) symmetry in the pseudospin limit. The conserved U(1) 


generator is given by [13] 
- (4, 0 
L, < ( 0 N . (9) 


where £, = U,£.U, and £, = (r x p);. In this case, the Dirac eigenstates of Hps are simultaneous eigenstates o of L,, 


with eigenvalue A, and of the total angular momentum generator J, = Š, + Lz, with eigenvalues Q = A + ñ —A £ i 
T ps A _ A RDS = 
L, O ALPI — A 7 aa” ; 
prs =A ez ps > 
J; $ A of) = OP i anu) 1 (10) 


Here 7) denotes additional quantum numbers that may be needed to specify the states uniquely. 

'The conventional method of labeling the eigenstates of axially deformed single-particle states in nuclei is to use 
the asymptotic quantum numbers (N,ns, A, Q), mentioned in the Introduction, that emerge in the limit of a non- 
relativistic axially-symmetric deformed harmonic oscillator with spin symmetry. For the relativistic axially-deformed 
harmonic oscillator with spin symmetry [29] the eigenfunctions can also be labeled by these quantum numbers. 
However, only the spatial amplitudes of the upper components of the doublet will necessarily have the nodes suggested 
by these quantum numbers, whereas the spatial amplitudes of the lower components may have different nodal structure. 
For spherically symmetric potentials a general theorem relates the nodal structure of the upper and lower Dirac 
amplitudes, and has been used to explain the non-relativistic radial quantum numbers characterizing pseudospin 
doublets in spherical nuclei [30]. A corresponding theorem for axially-deformed potentials in the pseudospin and spin 
limits of the Dirac Hamiltonian appears to hold under certain conditions which the relativistic harmonic oscillator 
satisfies, but which do not generally apply for realistic axially-symmetric potentials [31]. For the latter, only the 
quantum numbers A and Q in Eq. (10) are conserved in the pseudospin limit. The fact that the axial-symmetry of the 
potentials determines the ¢-dependence of the Dirac wave functions, leads to the following form for the relativistic 
pseudospin doublet eigenstates [14] 
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The two states in the doublet have the same pseudo-orbital angular momentum projection along the symmetry 


axis, A, but different total angular momentum projections Q = A — i and Q = A+ i. As seen from Eqs. (11a) and 
(11b), the pseudospin projection, ji = ti, and A are respectively the ordinary spin projection and ordinary orbital 
angular momentum projection of the non-vanishing lower component f; A (p, z). The corresponding dominant upper 


components ČRNA _i(p,z) and Ja 1(p,z) have orbital angular momentum projections A = A—1 and A’ = A+1 
$4497 3 7712 
respectively, hence A’ = A 4-2. Accordingly, Q = A + 1/2 and Q’ = A — 1/2 = ^ 4-3/2. These assignments agree with 
the non-relativistic pseudospin quantum numbers discussed in the Introduction. The generic label 7 in e x ü al?) 
replaces the harmonic oscillator labels N and ng, which are not conserved for realistic axially-deformed potentials in 
nuclei. 
In obtaining the expressions in Eq. (11) we have used the relations in Eq. (7), which for axially-deformed potentials 


read [14] 


faa a 02) = fa à 4052) 70, (12a) 
fixi = fux gb) = faa 2) > (12b) 
I7 i, (P z) = —95x (52) = 95A 0 2) » (12c) 
O À-1Y _ ð A-1 
ð 4 8 , A), 4 
Jz Iniz?) = (Z d 95 A ka 2) - (13b) 


We shall now test to see if the pseudospin symmetry conditions in Eqs. (12)-(13) are valid for realistic relativistic 
mean field eigenfunctions in deformed nuclei. 


IV. COMPARISON WITH REALISTIC RELATIVISTIC EIGENFUNCTIONS 


The single-particle energies and wave functions for !95Er are calculated by the relativistic Hartree theory with the 
parameter set NL3 in a Woods-Saxon basis [23, 24]. This method has been developed from relativistic theory in 
coordinate space [32, 33, 34] and has the advantage that it easily generalizes to include both deformation and pairing 
correction self-consistently. The pairing correlation is treated with the BCS approximation. These calculations lead 
to a theoretical average binding energy B/A = 8.107 MeV, a quadrupole deformation 6 = 0.3497 and a root mean 
square radius R = 5.376 fm, which reproduce the data well. For these realistic eigenfunctions the harmonic oscillator 
quantum numbers are not conserved, but the orbitals are labeled by the quantum numbers of the main spherical basis 
state in the expansion of the dominant upper component in the Dirac eigenfunctions. 

In Figure 1, the calculated single-neutron energies, € = E — M, for the pseudospin doublets in 'SEr are presented. 
From left to right, the panels correspond to the pseudo orbital angular momentum projection A = 1, 2, 3 and 4, 
respectively. The energy splitting between members of pseudospin doublets decreases as the single-particle binding 
energy, -e, decreases. For pseudospin doublets with binding energy larger than 5 MeV, the spin-up (pseudospin down) 
state is higher than the spin-down (pseudospin up) one. On the other hand, for the bound doublets with binding 
energy less than 5 MeV, the opposite is observed. 

Four pairs of neutron pseudospin partners are chosen to illustrate the relations given above. (i) The states [402]2 


and [404]2 (A = 3), which have a large energy splitting (about 2 MeV). The single-particle energies are respectively 
E4024 = —12.083 and E4047 = —14.160 MeV. (ii) The states [400]5 and [402]5 (A = 1), which have a small energy 


splitting (about 0.4 MeV). The single-particle energies are respectively € (499) i = —10.2073 and ziunas = —10.603 MeV. 


T 

z 
(iii) The states [501]2 and [503]2 (A — 2), which have a small energy splitting (less than 0.4 MeV). The single-particle 
energies are respectively zs = —1.349 and epggg3 = —0.9603 MeV. (iv) The states [510]2 and [512] (A = 1), 
which have a tiny energy splitting (less than 0.01 MeV). The single-particle energies are respectively £[510]3 = —3.8436 
and €[519)3 = —3.8378 MeV. 

Plots for the above four pairs of neutron pseudospin partners are shown in Figures 2, 4, 6, 8, as a function of z for 
three segments: p = 1, 3, 5 fm, and in Figures 3, 5, 7, 9, as a function of p for three segments: z = 1, 3, 5 fm. In each 
segment the top row displays the relationships between lower component amplitudes given in Eqs. (12a)-(12b), and 
the relationship between upper component amplitudes given in Eq. (12c). The bottom row displays the differential 
relationships between upper component amplitudes given in Eq. (13). 

From these figures, we can draw a number of conclusions. First, while the amplitudes [ Ki (p, z), m x (p, z) are 
not zero as predicted by Eq. (12a), they are much smaller than PEL z), PIC z). Furthermore, fia z) 
and fs x i002) have similar shapes as predicted by Eq. (12b). In Figures 2-3 for the [402]5/2, [404]7/2, A=3 
doublet there is some discrepancy in the shapes but the shapes become more equal as both the pseudo-orbital angular 
momentum projection A decreases (see Figures 4-5 for the [400]1/2, [402]3/2, A = 1 doublet) and the binding energy 
decreases (see Figures 6-7 for the [501]3/2, [503] 5/2, A = 2 doublet). 

The amplitude cÉ d (p,z) has the same shape as the amplitude 9x (p, z), in line with the prediction 


T 1 
of Eq. (12c), but they differ in magnitude. Again the discrepancy decreases as the pseudo-orbital angular 
momentum projection A decreases (compare Figures 2-3 for the [402]5/2, [404]7/2, A= 3 doublet with Figures 
8-9 for the [510]1/2, [512]3/2, A=1 doublet) and the binding energy decreases (compare Figures 2-3 for the 
[402]5/2, [404]7/2, A — 3 doublet with Figure 6-7 for the [501]3/2, [503]5/2, A — 2 doublet or Figures 8-9 for the 


[510]1/2, [512]3/2, A = 1 doublet). These amplitudes are much smaller than the other upper amplitudes, do mt 
ST 


The differential relation in Eq. (13a) between the dominant upper components, d. Al (p,z) and 9; i 
ALF ; 


p. 2), is 
AG (p, z) 


to the small upper components 9; A 41 (0,2). The shapes of the left-hand-side and of the right-hand-side of Eq. (13b) 
ALS 


are the same, but the corresponding amplitudes are quite different. Therefore, the differential relations in Eq. (13b) 
are less satisfied. These differences might partly originate from the differences in the magnitudes of the small upper 
components in Eq. (12c). 


1 
Dep 
well obeyed in all cases. The differential relations in Eq. (13b) relate the dominant upper components, g> 


[v 


V. SUMMARY 


We have reviewed the conditions that pseudospin symmetry places on the Dirac eigenfunctions. We have shown that 
the conditions on the lower amplitudes, Eqs. (12a)-(12b), are approximately satisfied for axially deformed nuclei. The 
differential relation between the dominant upper component amplitudes, Eq. (13a), is also approximately satisfied. 
However, both the relation between the amplitudes of the small upper components, Eq. (12c), and the differential 
equations, Eq. (13b), that relate the dominant upper components with the small upper components are not well 
satisfied. The pseudospin symmetry improves as the binding energy and pseudo-orbital angular momentum projection 
decrease, which is consistent with previous tests of pseudospin symmetry in spherical nuclei. 
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FIG. 2: Wave functions in (Fermi)~3/? as a function of z and p = 1, 3, 5 fm for the neutron pseudospin doublet [402]5/2 and 
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FIG. 3: Wave functions in (Fermi) ?/? as a function of p and z = 1, 3, 5 fm for the neutron pseudospin doublet [402]5/2 and 
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FIG. 4: As in Fig. 2 but for the neutron pseudospin doublet [400]1/2 and [402]3/2 (A =1) in 'SEr. 
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FIG. 5: As in Fig. 3 but for the neutron pseudospin doublet [400]1/2 and [402]3/2 (A =1) in 1° Er. 
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FIG. 6: As in Fig. 2 but for the neutron pseudospin doublet [501]3/2 and [503]5/2 (A =2) in 1° Er. 
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FIG. 7: Asin Fig. 3 but for the neutron pseudospin doublet [501]3/2 and [503]5/2 (A =2) in 'SEr. 
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FIG. 8: Asin Fig. 2 but for the neutron pseudospin doublet [510]1/2 and [512]3/2 (A =1) in 'SEr. 
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FIG. 9: As in Fig. 3 but for the neutron pseudospin doublet [510]1/2 and [512]3/2 (A =1) in 'SEr. 


